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Abstract 

In order to investigate the effect of inhomogeneities on the volume expansion of the universe, 
we study modified Swiss-Cheese universe model. Since this model is an exact solution of Einstein 
equations, we can get an insight into non- linear dynamics of inhomogeneous universe from it. We 
find that inhomogeneities make the volume expansion slower than that of the background Einstein- 
de Sitter universe when those can be regarded as small fluctuations in the background universe. 
This result is consistent with the previous studies based on the second order perturbation analysis. 
On the other hand, if the inhomogeneities can not be treated as small perturbations, the volume 
expansion of the universe depends on the type of fluctuations. Although the volume expansion 
rate approaches to the background value asymptotically, the volume itself can be finally arbitrarily 
smaller than the background one and can be larger than that of the background but there is an 
upper bound on it. 
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I. INTRODUCTION 



The standard Big Bang scenario is based on an assumption of the homogeneous and 
isotropic distribution of matter and radiation. This assumption then leads to the Robertson- 
Walker spacetime geometry and the Friedmann-Lemaitre (FL) universe model through the 
Einstein equations. This model has succeeded in explaining various important observa- 
tional facts: Hubble's expansion law, the content of light elements and the isotropic cosmic 
microwave background radiation (CMBR) [I]]. 

The CMBR conversely gives a strong observational basis for the assumption of homo- 
geneity and isotropy of our universe by its highly isotropic distribution together with the 
Copernican principle; we know that our universe was highly isotropic and homogeneous at 
least on the last scattering surface where CMBR comes from[Q. Hence, in the early stage 
of our universe, the linear perturbation analysis in the FL universe model is a powerful tool 
to investigate the dynamical evolution of our universe [0]. 

In order to perform the perturbation analysis, we need an appropriate background uni- 
verse model, i.e., information of the Hubble parameter, the density parameter and further 
the equation of state of the matter, radiation and so on, in the real universe. To fix the 
background universe, we use the observational data in the neighborhood of our galaxy. Espe- 
cially the Hubble parameter is determined from the data about the distance-redshift relation 
within lOO/i^Mpc except for the type la supernova^, E|. However, the universe in a region 
within 100/i _1 Mpc is highly inhomogeneous and hence there are non-trivial prescriptions to 
identify the present inhomogeneous universe to the homogeneous and isotropic FL universe 
model. If those procedure are not appropriate, we might miss finding the correct background 
universe. 

It is often stated that the spatially averaged observational data in the vicinity of our 
galaxy are recognized as those of the background FL universe model. The Hubble parameter 
determined by the observed distance-redshift relation in our universe is regarded as the 
expansion rate of the volume of the region co-moving to matter. There are several researches 
for the effects of inhomogeneities on the volume expansion of the universe]^, 0, H, ^, [It], [TT 



12], O, [14], [I5(|. Especially, Nambu applied the renormalization group method to the second 



order cosmological perturbation theory and claimed that the expansion of the dust filled 



universe is decelerated by the inhomogeneities fl2" 
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In the real universe, this back reaction effect might be very small ||. However in order 
to get deeper insight into the dynamics of the inhomogeneous universe, we will consider the 
situation in which the back reaction of the inhomogeneities seems to be effective. For this 
purpose, we consider the Swiss-Cheese universe model and investigate the volume expansion 
in it. The original Swiss-Cheese universe model is constructed by choosing non-overlapping 
spherical regions in the background homogeneous and isotropic dust filled universe and then 
replacing these regions by the Schwarzschild space-time whose mass parameter is identical 
with the "gravitational" mass of the dust fluid in the removed region. On the other hand, 
in this article, we consider a modified version; we first remove spherical regions from the 
homogeneous and isotropic dust filled universe and then fill these regions with spherically 
symmetric but inhomogeneous dust balls. A spherically symmetric inhomogeneous dust ball 
is described by the Lemaitre-Tolman-Bondi (LTB) solution which is an exact solution of the 
Einstein equations, and hence by this procedure, we obtain an exact solution of Einstein 
equations, which represents an inhomogeneous universe. Using this solution, we can study 
non-linear effects of inhomogeneities on the volume expansion of the universe without use 
of perturbation analysis. 

In the LTB solution, shell crossing singularities are generic. Since the LTB solution is no 
longer valid after the occurrence of the shell crossing, we need to change the treatment if 
it occurs. As a crude approximation to describe the dynamics after the shell crossing, we 
adopt a model in which the shell crossing region is replaced by a spherical dust shell. 

This article is organized as follows. In section II, we explain how to construct modified 
Swiss-Cheese universe models which are studied in this article. We investigate the volume 



expansion rate in the case of small perturbations in section III and in the highly inhomoge- 



neous case in [lv]. In section [lv], an alternative model is also constructed which describes the 
universe after the shell crossing and investigate the dynamics of this model. Finally, section 
[V] is devoted to summary and discussion. 

We use the units in which c — G—1 throughout the paper. 



II. MODIFIED SWISS CHEESE UNIVERSE MODEL 

In this section, we give a prescription to construct MSC universe model. First we con- 
sider a Einstein-de Sitter (EdS) universe and then remove spherical regions from it; these 
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removed regions should not overlap with each other. Next these regions filled with inho- 
mogeneous dust balls with the same radii and the same gravitational mass as those of the 
removed homogeneous dust balls. In this MSC universe model, each inhomogeneous region 
is described by the Lemaitre-Tolman-Bondi (LTB) solution which is an exact solution of 
Einstein equations. 

LTB solution describes the dust filled spherically symmetric spacetime. Adopting syn- 
chronous and co-moving coordinate system, the line element is written as 

ds 2 = — dt 2 + r )ijdx l dx^ 

= -dt 2 + ^&^- d X 2 + Y 2 (t, X )(d6 2 + sin 2 9dip 2 ), (1) 

where the prime ' denotes the differentiation with respect to the radial coordinate x- in this 
coordinate system, components of 4-velocity u a of a dust fluid element are 

u a = (1, 0, 0, 0). (2) 

The stress-energy tensor T ab is then given by 

T ab = p(t,x)6 a 6° b , (3) 

where p(t, x) is the rest mass density of the dust. 

Einstein equations lead to the equations for the areal radius Y(t, x) and the rest mass 
density pit, x) °f t ne dust; 

+ (4) 

where k(x) and M(x) are arbitrary functions and the dot ' denotes the differentiation with 
respect to t. 

We set M(x) as 

Mix) = ^X 3 , (6) 

where po is a non-negative arbitrary constant. The above choice of M(x) does not loose any 
generality. Eqs. are invariant for the rescaling of the radial coordinate x-, 

X^X = X(x)- (7) 
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Considering this property, we can choose above form of M(%) as long as pY' > 0. 

The solutions of eq. (Q) are given as follows: 
In the region where k(x) > 0, 

4vrp 0l 



Y 
to(x) 



3k 
'3k 3 / 2 



[1 - COST])x, 

(r) — sin rj); 



in the region where k(x) = 0, 



in the region where k(x) < 0, 



Y 



67ip {t - t (x)Y 



1/3 



XI 



\ = _ , ^° (cosh — 



t-to(x) 



3\k\ 
4vrp 



(sinhr] — r/), 



(8) 
(9) 

(10) 

(11) 
(12) 



3|fc| 3 / 2 

where to(x) is an arbitrary function. Note that to(x) is the time when a shell focusing 
singularity appears, where 'shell focusing singularity' means Y = for x > and V = at 
X = 0. In this article, we consider a region of t > £q; an d hence the time t = t corresponds 
to the Big Bang singularity. Here, we focus on a case of to — 0, i.e., simultaneous Big Bang. 
For simplicity, we consider the simplest version of MSC universe models shown in fig. [I]; 
there is only one inhomogeneous spherical region at the center in each identical cubic region 
fl We focus on only one cubic co-moving region Q (fig. 0). In this article, we consider a 
following model. Assuming < xi < X2 < X3 < Xsc, 



k(x) 



k 
kg 

h 

h 
[2x 2 



(* 2 ~ ^) 2 2 2 

2 2 ^ Al ^ A2 

XI - X2 

(xi + xl) 



(xi + xl) 




1 



for < x < Xi 

for xi < X < X2 

for X2 < X < X3 

for X3 < X < Xsc 



(13) 



where ko is constant. In order to guarantee 1 — x 2 k > 0, the following inequality should 
hold, 



«:=Y(X? + X!)<1. 



(14) 
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FIG. 1: The MSC universe model. Each shaded region represents the inhomogeneity and is de- 
scribed by LTB solution. 




FIG. 2: One cubic region 0, of the MSC universe model. I and Xsc are the co-moving scales of this 
cubic region and the inhomogeneous region respectively. 
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We consider two cases; one is the case of k > and the other is the case of k < 0. Then 
we investigate the volume expansion rate of a cubic co-moving region Q. The volume V is 
defined by 



V(t) := / yffdPx, (15) 
Jn 

where 7 is the determinant of the spatial metric 7^. The volume expansion rate is defined 
as V/V. 



III. THE CASE OF SMALL FLUCTUATIONS 



In a region where 



< ' 1 < 1 
2vrp 



is satisfied, the areal radius Y(t, x) is written in the form of power series as 

1 3 



Y(t, x) = a(t) X ( 1 - j^e - ^e 2 ) + 0(e 3 ), 



where 



, 9t V /0 , ( !2i Y 
e(t,x):=[- k= hTTTT X% 



and 



2 / 3 / Yit \ 2 ^ 3 
a(t) := (6np t 2 ) 1/3 . 



(16) 



(17) 



(18) 



Further, we consider the case of \k\x 2 <C 1, which has been studied by Nambu by the second 
order perturbation analysis. The components of the metric tensor are written as 



9xx = a 



( M\ 



2/3 



1 d 



— l e-— ^(xe)+0(e 2 ) 
\l2t) lOrfx 



2 2 
a X 



2 2 • A n 

g<p<p — a X sm v 



i-^ + o( e 2 ) 



1-— e + 0(e 2 ) 
10 V ; 



(19) 
(20) 
(21) 



From the above equations, it is easily seen that in the limit e — > with t fixed, the metric 
tensor becomes that of EdS. Since the outside region is EdS universe, e should vanish at the 
boundary x — Xsc by the continuity of the metric tensor. 



To compare our result with the study by Nambu|T2||, we impose conditions that spatial 
averages of the Cartesian components of the metric tensor and of the density agree with 
those of the background EdS universe up to the first order of e. 

The spatial average of a quantity F is defined as follows: 



(F) ,- 



(22) 



d 3 x) / Fd 3 x. 

We consider a Cartesian coordinate system (x, y, z) which is related to the spherical polar 
coordinate system (x, 9, if) in the ordinary manner as 

x = x sin 9 cos <p, 
y = x sin 9 sin ip, 
z = x cos 9. 

Then the spatial averages of Cartesian components of the spatial metric 7^ are obtained as 



(jij) = a 2 



' 4tt [*° \ ( M\ 2/3 2 1 d . 3 , . , 



S lJ + 0(e 2 



(67rPot2)2/3 I 1 + £ kx4dx 1 Sij + 0{e ^ 



(23) 



where we have used e = at x — Xsc- Therefore we should define the scale factor a B (t) of 
the background EdS universe as 



47T f Xsc N 1/2 

a B {t) := a(t) [l + —J^ k X 4 d X 



(24) 



The above equation means that although the outside homogeneous region has the same 
geometry of EdS universe, it does not agree with the background EdS universe as long as 

"Xsc 



kx"d X + 0. 



(25) 



The rest mass density p of the dust is written in the form of a power series with respect 
to e as 



6nt 2 { 20x 2 d X 
The spatial average of p is obtained as 

1 



(p) 



6vrt 2 



+ 0(e 2 ). 



(26) 



(27) 



Hence the background energy density p B is defined by 

1 



Pb(*) := 



6irt 2 ' 

Here note that the background Hubble equation 

2 



(28) 



>7T 



-P 



B) 



(29) 



holds. 



By eq. (|16|), the 3-dimensional volume element ^7 is written as 



y/l = a 3 * 2 sin 



+ 0(e 3 ). (30) 



Using the above equation, the volume V defined by eq. ( jT5| ) is obtained as 



where 



y(t) = a 3(t) (£ 3 + V 1 + K 2 t 2/3 ) + 0(e 3 ), 



(31) 



Vo 



'o 3£ 6 \Jo 

20 \ 2^p~ ) 



k\ A d X < 0. 



(32) 
(33) 



In eq. fl3~T|), a|£ 3 is the 3-dimensional volume measured by background EdS geometry and ex- 
tra terms come form inhomogeneities. These terms do not include first order perturbations, 
but come from second order perturbations. 
The volume expansion rate is given by 



(34) 



The first term of the R.H.S. in the above equation corresponds to the background part. On 
the other hand, the second term implies that the back-reaction of inhomogeneities decelerates 
the volume expansion. It is worthwhile to note that this result does not depend on the 
detailed functional form of fc(x). 
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IV. THE CASE OF NON-LINEAR FLUCTUATIONS 

In this section, we study the cases where \k\x 2 is not necessarily much smaller than 
unity. As in the case treated in the previous section, in order to specify an effect due 
to inhomogeneity, we need a background homogeneous cubic region to be compared with 
an inhomogeneous one. However the background homogeneous universe introduced in the 
previous section is not appropriate for the non-linear case; for example, too small k makes the 
background scale factor a B defined in eq. (J23|) negative. In order to introduce an appropriate 
background, we consider the rest mass Mr defined by 



where g is the determinant of the metric tensor of spacetime. Mr is a conserved quantity 
by virtue of the continuity of the rest mass density, d a (pu a y/—g) = 0, where d a is a partial 
derivative. We introduce a background as a cubic region with the same rest mass as the 
corresponding inhomogeneous cubic region. Note that in general, the rest mass of the dust 
Mr in a cubic co- moving region disagrees with that of the original EdS universe. Hence a 
cubic region of the original EdS universe is not background. 

The rest mass density p B and scale factor a B of the background are introduced in the 
following manner, 



Even if we specify Mr, p B and a B are not fixed completely; we need one more condition. 
Here we impose a condition in which the volume V(t) approaches to the volume a 3 (t)£ 3 of 
the background cubic region for t — > 0. In the limit of t — > 0, the volume V behaves as 




(35) 



p B a 3 J 3 : = Mr. 



(36) 




(37) 
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Hence the background scale factor a B is given by 



1/3 



a B (t) := a { 1 + % [*" [ j^== ~ 1 ] X*d X } ■ (38) 



o 



Here note that in the case of x 2 \k\ ^ 1> the above definition of a B agrees with eq. ([24|) up 
to the first order of x 2 k- From eqs. flSHj), (|36|) and (|38"D, we find that the background rest 
mass density p B is completely the same as eq. (p8|). We can easily see that the background 
Hubble equation (p9|) also holds. 

From eq. ([|), we can easily see that if Y' vanishes, the rest mass density p becomes infinite 
and hence a singularity forms there. This is called shell crossing singularity. Hellaby and 
Lake showed that a necessary and sufficient condition for the appearance of a shell crossing 
singularity is |IB | 



k' > for the region k > 0, (39) 
(x 2 k)' > for the region k < 0. (40) 

In the case of ko > 0, the first condition does not hold and hence a shell crossing singularity 
does not appear. On the other hand, in the case of ko < 0, shell crossing singularities Y' — 
necessarily appear since (x 2 ^)' > in the region xs < X < Xsc- 



A. The Case of k > 

To estimate the volume V, we rewrite it in the form, 

^° 3 ( <3 -H +4 *f vSf (41) 

where Y sc := a(t)xsc- From eq. @, we can see that K vanishes at = 2tt by gravitational 
collapse, and hence by substituting i] = 2n into eq. (§), the singularity formation time 
^ = ^sg(x) is obtained as 

Denoting the inverse function of t 8S (x) by Xs g(t), the region of < x < X sg(t) < Xsc has 
already collapsed at time t larger than 8ir 2 p /3k^ 2 . Xsg(t) approaches to Xsc for t —>■ oo 
asymptotically. Here note that in the integrand of eq. (0), x — Xag(t) at Y = and X — Xsc 
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at Y = Y BC . Since x 2 k(x) is decreasing function with respect to x in the region xz < X < Xsc 
and vanishes just at x — Xsc, we find that < x 2 k(x) < XsgMXsg) holds in the integrand 
at sufficiently large t. This means that x 2 k{x) also approaches to zero asymptotically, since 
XsgMXsg) -> Xsc^(Xsc) = for t -> oo, and thus 

yxL)+4vr^ yS VW = a 3 £ 3 . (43) 

This equation means that the volume expansion rate approaches to the background value 
asymptotically, i.e., 

V a d R . 

- — > 3- = 3— for t -> oo. 44 
1/ a a B 

However the volume itself may approach to much different value from the background one 
a 3 £ 3 . \J\ — x 2 k can be made arbitrarily smaller than unity in the region xi ^ X < X3j m 
the limit of k — > 1, a/1 — x 2 /c — > in this region (see eqs. (|13|) and (0)). Thus, if we set k 
to be very close to unity, we obtain 

1/3 

In this case, the volume V approaches to the value much smaller than the background one, 
asymptotically. The volume expansion is also much different from that of the background 
in the intermediate stage (see fig. ||). 



. 4tt f Xsc X 2 dx i 
«B~«k 7=== » a - ( 45 ) 



B. The Case of k < 

As mentioned in the above, shell crossing singularities Y' = necessarily appear in this 
model. Before it, the volume expansion is shown in fig. For k(x) < 0) background scale 
factor a B cannot differ from original one very much. We plot the ratio of the volume V(t) 
to a 3 £ 3 and a^i 3 . We find that the inhomogeneities decelerate the volume expansion before 
shell crossing singularity appears. 

Here, we investigate the volume expansion rate after the appearance of this shell crossing 
singularity. The structure formed by the shell crossing depends on what is approximated 
by the dust matter. In case that the dust matter is extremely cold fluid, a shock wave will 
form after the shell crossing. If the dust matter consists of collisionless particles, a spherical 
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; EdS background 




1 2 

Time 



FIG. 3: Volume expansion with positive k(x)- The dotted line is the temporal evolution of the 
co-moving volume measured by the background EdS geometry and the dashed line is that measured 
by the outer EdS geometry. The time is set to unity, when the scale of the spherical inhomogeneous 
region a(t)xsc agrees with that of the horizon scale a(t)/a(t). The scale factor at t = 1 is also set 
to unity and the co- moving scale £ is set to 4x sc - Hence the co- moving volume measured by the 
outer EdS geometry at t = 1 is £ 3 (= 6 3 = 216). 

wall will form. When the width of the shock wave or of the wall is much smaller than the 
radius of it, we can treat the shock front or the wall as a timelike singular hypersurface, 
where 'timelike' means that the unit normal vector n a to the hypersurface is spacelike, i.e., 
n a n a = 1. 

A timelike singular hypersurface is characterized by its surface-stress-energy tensor de- 
fined by 

S ab := lim / T cd h c a h d b dx (46) 

J-e 

where x is a Gaussian coordinate (x = on the hypersurface) in the direction of the normal 
vector n a , and h c a := 5 c a — n a n c is a projection operator. A timelike singular hypersurface 



13 



1.02 " 



o 



1 



Vmsc/Vb 



0.98 - 



Vmsc/^out 



J_L 







0.03 



0.06 



0.09 



Time 



FIG. 4: Temporal evolution of the ratios before shell crossing singularity appears. The time is 
set to unity, when the scale of the spherical inhomogeneous region a(t)xsc agrees with that of the 
horizon scale a(t)/a(t). Dashed line corresponds to the shell crossing time. V MSC is the volume of 
the co-moving cubic region Q in modified Swiss Cheese universe. V B is the volume of O measured 
by the background EdS geometry, i.e., V B = a|^ 3 . V OVT is the volume of Q measured by the outer 
EdS geometry, i.e., V^ UT = a 3 ^ 3 . 

with a surface-stress-energy tensor of the form 



is called a world sheet generated by a trajectory of a dust shell, where a is the surface-energy 
density and v a is the 4- velocity of an infinitesimal surface element of the dust shell. Hereafter 
we focus on this case. 

In order to get an insight into the dynamics after shell crossing, we study the volume 
expansion of a cubic region with a spherically symmetric dust shell. In the case of ko < 0, 
k(x) is negative in the LTB region; (x 2 k)' < in the inner region < x ^ X3> while 
{x 2 k)' > in the outer region X3 < X < Xsc- I n accordance with eq. ([5]), shell crossing 
necessarily occur in the outer region and we assume that this region collapses into a dust 



Sab = crv a v b , 



(47) 
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shell. Hence we focus on a situation in which (x 2 k)' < holds inside the dust shell. The 
model is constructed by enclosing the interior LTB region by a spherically symmetric timelike 
singular hypersurface (see fig. |]). 

The analysis by using the dust shell model works only before the dust shell reaches the 
boundary of the cubic co-moving region Q. When the dust shell reaches the boundary, it 
collides with other dust shells centered in surrounding cubic regions. In this article, we do 
not consider the dynamics after collisions of the dust shells and assume that there is enough 
time before collisions of dust shells. 

The interior region of the dust shell is described by the LTB solution with a line element, 

ds 2 _ = -dt 2 _ + ^'if^ dxt + Y 2 (*_ , x- ) (dO 2 + sin 2 6d<p 2 ) , (48) 

where the prime ' means the derivative with respect to X- ■ The equation for the areal radius 
Y is given by the same equation as eq. (f|), i.e., 

= + ?^fc>, (49) 

where the dot ' means the derivative with respect to £_. Here we focus on the late time 
behavior of the volume expansion. In this case, since Y is monotonically increasing with re- 
spect to t-, the "gravitational potential" 2M(x-)/Y becomes much smaller than — xtk(x-) 
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and hence we ignore this potential term in eq. (f49|) . This approximation corresponds to an 
assumption in which the interior region of the dust shell is described by Minkowski geometry. 
Then the solution for Y is easily obtained as 

r(t_,x-) = {-x^( X -)} 1/2 t-. (so) 

The line element of the outer region is written as 

ds\ = -dt\ + a 2 (t+) {dx 2 + + xl {d0 2 + sin 2 0dp> 2 ) } , (51) 
and Einstein equations reduce to the Hubble equation as 

{ f = AMx ± ) j m 

where M + (x+) is related to the rest mass density p + of the outer region as 

M + ( X+ ) = —p + a 3 xl- (53) 

Note that by virtue of the spherical symmetry, angular coordinates, 9 and (p, are common 
for both interior and exterior regions of the dust shell. 

A physically and geometrically clear prescription to treat a timelike singular hypersurface 
has been presented by Israel||17||. In his prescription, junction conditions on metric tensor 
across the singular hypersurface lead to equations to determine the singular hypersurface 
itself, i.e., the equation of motion of a dust shell in our case. Using his prescription, dynamics 
of a vacuum void surrounded by a spherical dust shell in expanding universe has been 
analyzed by Maeda and SatoHl8|. We can use their results since the situation considered 
here is completely the same as theirs. 

By virtue of the spherical symmetry of the system considered here, the trajectory of the 
dust shell is given by 

£_ = t s _(t + ), x± — Xs±(t + ), = constant and (p = constant, (54) 

where the time coordinate t + in the exterior region has been adopted as an independent 
temporal variable. The areal radius R of the dust shell is then given by 

R(t+) := a (t+) Xs+ = Y- (*-, X*-) (55) 
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Maeda and Sato derived a differential equation for the areal radius R of the dust shell 
= ^ {- (1 + ^ + 2V 2 + V 2 ) + (1 - 4V 2 )(1 + ZVVh + ^l) 1 ^} , (56) 



where 



if := 



a(t. 



Vh := Fit!, 



(57) 

(58) 
(59) 



Using a solution of eq. flSED, the coordinate radius Xs+ of the dust shell in the outer region 
is given by 

R 



Xs+ 



a(t+)' 



(60) 



Equations for t s _ and Xs- on the dust shell are given by 



dt* 



{l-xlk( X s-)} 1/2 ^l-a 2 (t + ) 

- {-xik( Xs -)} 1/2 



dt A 



dXs 
dt A 



dR 



1/2 



dXs 



1 



dR ■ , . (it s 

Y-{t s -,Xs-)- 



(61) 
(62) 



dt + YL{ts-,Xe-) \dt+ V " B ' dt A 
where the dot ' and the prime ' are derivatives with respect to £_ and x_, respectively. The 
volume V\ n inside the dust shell is written as 

r-xs- v'v 2 r R y 2 



Vin = 4tt 



2ntt 



d\ = Ait / ===== 

o V 1 - x 2 Hx) Jo V 1 + Y2 /t't 



■dY 



R 



R 



In 



R 



R 



(63) 



where we have used eq. ( BTiD in the second equality to estimate x 2 k(x) i n the integrand. We 
consider a normalized volume V{ n by the volume of the removed homogeneous dust ball in 
the original EdS universe, 



3 (U 



AnR 3 /3 2\R 



R 



R 



R 



(64) 
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Vin is monotonically decreasing function of -R/t s _; it approaches to unity in the limit of 
R/t s - —>■ 0, while it vanishes in the limit R/t s - — > oo (see fig. |]). In order to see temporal 
behavior of Vi n , we need to solve eqs. ( [56"D and (|6T|). 

1.0 f ' ' - 




40 80 



R/t s - 

FIG. 6: Behavior of V in function of R/t B —. 



For sufficiently large £+, the areal radius i? of the dust shell becomes much smaller than 
the cosmological horizon scale H~ x . In this case, the the motion of the dust shell is well 
described by the Newtonian approximation. Maeda and Sato showed that for sufficiently 
large t + , R behaves as|T8| 

R(t+) OC t + (!5+VT7)/24 „ t+ 0.797 „ 4 _0.797 (g5) 

Hence we find that R/t B - oc t s _ -0 ' 203 — > for t s _ — > oo. Using this result and eq. we 
find Vin — > 1 for £ s _ — > oo, and hence for t + — ► oo, 

V — ► a 3 (i+)£ 3 . (66) 

This equation means that the volume expansion rate approaches to the background value 
asymptotically, i.e., 

^ — . 3^ = 3^ for t - oo. (67) 
V a a n 
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The effect of inhomogeneities on the volume expansion rate vanishes after the dust shell 
becomes much smaller than the horizon scale H^ 1 (see fig. [7p. However, as in the case of 
ko > 0, the volume itself is different from the background value a 3 £ 3 . By eq. (p8|), we find 
that the asymptotic value of V is larger than the background value a 3 £ 3 . However it should 



VmSc/Vb 




0.16 1 2 

Time 

FIG. 7: Temporal evolution of the ratios after shell crossing singularity appears. The time is set 
to unity, when the scale of spherical inhomogeneous region a(t)x S c agrees with that of the horizon 
scale a(t)/a(t). V MSC is the volume of the co-moving cubic region Q in modified Swiss Cheese 
universe. V B is the volume of measured by the background EdS geometry, i.e., V B = a 3 £ 3 . V 0VT 
is the volume of Q measured by the outer EdS geometry, i.e., V 0U t = a 3 £ 3 . 

be noted that there is an upper limit on the asymptotic value of V. Since k can be arbitrarily 
small in this model, a/1 — x 2 k can t> e niade arbitrarily larger than unity except at x = 
and x = Xsc- Hence we obtain 

/ 47]- X 1 / 3 / 7T\l/3 

a B (*+) > a(t + ) 1^1 - — xij > a{t+) (1 - -J , (68) 
where the last inequality is obtained by setting x sc = £/2. Using this inequality, we obtain 
a 3 B {t + )£ 3 < a 3 {t + )£ 3 < (l - ^j' 1 a 3 B {t + )£ 3 ~ 2.10 x a 3 B {t + )£ 3 . (69) 
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In contrast with the case of ko > 0, the volume itself can not be so different from the 
background value. 

V. SUMMARY AND DISCUSSION 

We have investigated an effect of inhomogeneities on the volume expansion in modified 
Swiss-Cheese universe model. We considered two cases; the inhomogeneities collapse into 
black holes (ko > 0), while inhomogeneities expands faster than the background volume 
expansion (k < 0). When inhomogeneities can be treated as perturbations of Einstein-de 
Sitter universe, the volume expansion is decelerated due to the second order contribution of 
the perturbations in both models. This result agrees with Nambu's second order perturba- 
tion analysis. 

Although the choice of background homogeneous universe is straightforward in the case 
of \k\x 2 <C 1, it is not in the case of non-linear situation. We introduced the background 
homogeneous universe in order to satisfy the conditions; the cubic region in the background 
homogeneous universe with the same rest mass have the same evolution of the volume as 
that of corresponding region in the modified Swiss-Cheese universe in the limit of t — ► 0. 

In the case of non- linear fluctuation with ko > 0, we find that the volume expansion 
rate approaches to that of background universe asymptotically. Since the modified Swiss- 
Cheese model is an exact solution of Einstein equation, we can obtain the precise behavior 
of the volume expansion. We set it so that for t — > 0, the temporal evolution agrees with 
background one. Then we found that for t — > oo, the temporal evolution agrees with that of 
outer EdS universe. We can see these asymptotic behaviors analytically but we have to rely 
only on the numerical method to obtain the behavior of the intermediate stages. From the 
result of the numerical calculation (fig. |3[), we found that volume expansion is decelerated 
by the inhomogeneities. This behavior coincides with a result obtained by Nambu. We 
note that his result is based on the perturbation theory but our result is not. In our highly 
non-linear example, the volume expansion rate becomes negative at the intermediate stage. 
This result may be the characteristics of the effect of non-linear fluctuations which can not 
be treated by the perturbation method. 

In the case of ko < 0, the shell crossing singularity appears in the inhomogeneous region. 
In this article, we assumed that a spherical dust shell forms after the shell crossing singu- 
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larity appears. We find that the inhomogeneities decelerate the volume expansion before 
it appears (fig. This is consistent with the previous results by Nambu but the inhomo- 
geneities cannot be treated as the perturbation of homogeneous universe. After a spherical 
dust shell forms, the volume expansion rate approaches to that of the background universe 
asymptotically. Our dust shell universe model (fig. |5|) is a crude approximation. Therefore 
the behavior of the volume expansion obtained (fig. |7]), especially at the early stage after 
the shell crossing, contains the influences of this crudeness. But the asymptotic behavior 
(t — > oo) might be free from this approximation. 

We fixed the form of k(x)- Temporal evolution of the volume depends on the form of 
k(x) but the asymptotic behavior may be independent of the form of fc(x); 

V — >a 3 £ 3 for t -> oo. 
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